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ABSTRACT 


Thermostructural analysis was performed on a heated titanium honeycomb-core sandwich panel. The 
sandwich panel was supported at its four edges with spar-like substructures that acted as heat sinks, 
which are generally not considered in the classical analysis. One side of the panel was heated to high tem- 
perature to simulate aerodynamic heating during hypersonic flight. Two types of surface heating were 
considered: (1) flat-temperature profile, which ignores the effect of edge heat sinks, and (2) dome- 
shaped-temperature profile, which approximates the actual surface temperature distribution associated 
with the existence of edge heat sinks. The finite-element method was used to calculate the deformation 
field and thermal stress distributions in the face sheets and core of the sandwich panel. The detailed ther- 
mal stress distributions in the sandwich panel are presented, and critical stress regions are identified. The 
study shows how the magnitudes of those critical stresses and their locations change with different heat- 
ing and edge conditions. This technical report presents comprehensive, three-dimensional graphical dis- 
plays of thermal stress distributions in every part of a titanium honeycomb-core sandwich panel subjected 
to hypersonic heating on one side. The plots offer quick visualization of the structural response of the 
panel and are very useful for hot structures designers to identify the critical stress regions. 

NOMENCLATURE 



m 

SPAR 

S81 

T 


length of sandwich panel, in. 

width of sandwich panel, in. 

modulus of elasticity of sandwich face sheets, lb/in 2 

effective modulus of elasticity of sandwich core in x-direction, lb/in 2 

effective modulus of elasticity of sandwich core in y-direction, lb/in 2 

effective modulus of elasticity of sandwich core in z-direction, lb/in 2 

beam element for which the intrinsic stiffness matrix is given 

quadrilateral combined membrane and bending element 

effective shear modulus of sandwich core in xy-plane, lb/in 2 

effective shear modulus of sandwich core in xz-plane, lb/in 2 

effective shear modulus of sandwich core in yz-plane, lb/in 2 

depth of sandwich panel, in. 

depth of sandwich core, in. 

joint location (or grid point or node) of finite-element model 
number of deformation half waves in ^-direction 
structural performance and resizing finite-element computer program 
hexahedron (or brick) element 
temperature, °F 

temperature of lower face sheet, °F 



T u temperature of upper face sheet for the flat-temperature profile, or temperature of 

upper-face-sheet plateau zone for the dome-temperature profile, °F 

t s thickness of sandwich face sheets, in. 

w deflection at arbitrary point of middle plane of sandwich panel, in. 

w max maximum deflection at center of sandwich panel, in. 

x, y, z rectangular Cartesian coordinates, in. 

x' shifted x-coordinate (x' = x + a/2), in. 

a coefficient of thermal expansion of solid plate or sandwich face sheets, in/in-°F 

« CJC coefficient of thermal expansion of sandwich core in x-direction, in/in-°F 

c i coefficient of thermal expansion of sandwich core in y-direction, in/in-°F 

ct C2 coefficient of thermal expansion of sandwich core in z-direction, in/in-°F 

oc l coefficient of thermal expansion of sandwich face sheets at temperature T in/in-°F 

cl u coefficient of thermal expansion of sandwich face sheets at temperature T , in/in-°F 

AT temperature differential between upper and lower face sheets (AT -T - T t ), °F 

y xz , V yz transverse shear strains of sandwich panel in xz- and yz-plane, in/in. 

v Poisson ratio of sandwich face sheets 

v cxy ’ v cxz ’ v cyz P°i sson ratios of sandwich core 

p hc density of sandwich core, lb/in 3 

o x normal stress in x-direction, lb/in 2 

o y normal stress in y-direction, lb/in 2 

x shear stress in Ay-plane, lb/in 2 

T xz ’ T yz transverse shear stresses in sandwich core in xz- and yz-planes, respectively, lb/in 2 


INTRODUCTION 

A sandwich panel fabricated with titanium face sheets bonded to titanium honeycomb core through 
enhanced diffusion bonding process is a potential candidate for application to hypersonic aircraft outer 
skin structural panels (ref. 1). This type of sandwich structure can operate at elevated temperature levels 
approaching 1000 °F. When applied as a structural component of hypersonic flight vehicles, this type of 
sandwich panel is fastened to relatively cool substructures that act as heat sinks. Even under uniform 
surface heating, the induced panel surface temperature distribution could be nonuniform because of those 
edge heat sinks. Most analyses do not include the heat sink effects because of added mathematical com- 
plexity (refs. 2 through 8). The heated sandwich surface temperature profile is generally a truncated dome 
shape, with temperature nearly constant in the central plateau zone, tapering down toward the cooler 
edges. Because the panel is supported by relatively cool substructures and constrained from free 
expansion, considerable thermal stresses could build up in the panel. The most critical stresses are the 
compressive stresses. Excessive magnitude of compressive stress built up in the heated face sheet could 
cause thermal bending caused by thermal moments; thermal buckling; thermal yielding; thermal creep; 
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thermal crack after cooling down; and other effects. One-sided heating, under certain temperature 
profiles and edge conditions, also could induce high-intensity transverse shear stress in the sandwich core 
near the panel comer, which could cause potential shear debonding between the face sheets and the sand- 
wich core. Thus, loss of structural integrity could result. 

Ko and Jackson conducted extensive studies, in recent years, concerning the mechanical and thermal 
buckling characteristics of titanium sandwich panels (refs. 2 through 6) and metal-matrix composite 
sandwich panels (refs. 7 and 8). Extensive information about thermomechanical buckling characteristics 
of such sandwich structures have been documented (refs. 2 through 8). To fully understand the thermo- 
structural response of the sandwich panels in actual applications under which the panel is constrained by 
the substructures and subjected to one-sided heating, detailed thermal stress analyses are needed to iden- 
tify the critical stress regions. 

This report presents the results of finite-element thermal stress analyses of the sandwich panel under 
different one-sided heating conditions and edge constraints. The detailed deformation fields and thermal 
stress fields generated in the sandwich panel are presented graphically for easy visualization of the 
critical stress regions. 


DESCRIPTION OF PROBLEMS 

Figure 1 shows the honeycomb-core sandwich panel, which has length a , width b, and depth h (depth 
of sandwich core h Q ). The upper and the lower face sheets have the same thickness of t . The panel is 
subjected to one-sided heating of A T-T — the temperature differential between the upper-face-sheet 
temperature T and the lower-face-sheet temperature T ^ The temperature differential AT has two types 
of profiles (i.e., distributions): flat (fig. 2), for which AT is constant over the panel surface, and dome 
shaped (fig. 3), for which AT is constant only in the panel central region and decreases linearly to zero at 
the panel edges. The flat temperature profile heating is for the case when the heat sinks at the panel edges 
are neglected. The dome-shaped temperature profile heating is for the case when there exist cooler sub- 
structures at the panel edges. The temperature profile in figure 3 is actually a truncated pyramid and 
approximates the actual dome-shaped temperature profile (fig. 4), measured during a thermal ground test 
of a titanium sandwich panel heated on one side at 10 °F/sec heating rate (ref. 9). 

In the thermostructural analysis, the extensional and bending stiffnesses of the sandwich panel were 
provided by the two face sheets, and the transverse shear stiffness by the sandwich core. The sandwich 
panel was supported under four edge conditions to study different thermal deformation and thermal stress 
fields generated in the panel: 

1. 4S fixed-edge condition — conventional simply supported edge condition in which the four edges 
cannot move in the x-, y-, or z-directions. 

2. 4S free-edge condition — simply supported edges in which the four edges can move freely in the x- 
and y- directions only. 

3. 4C fixed-edge condition — conventional clamped edge condition in which the four edges have zero 
slopes and cannot move in the x-, y-, or z-directions. 

4. 4C free-edge condition — clamped edges with zero edge slopes in which the four edges can move 
freely in the x- and y-directions only. 
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FINITE-ELEMENT ANALYSIS 


This section describes the finite-element models and numerical input values used in the analysis. 

Finite-Element Modeling 

The structural performance and resizing (SPAR) finite-element computer program (ref. 10) was used 
in the thermostructural analysis of the sandwich panel. Because the panel is symmetrical with respect to 
the x- and y-axes (fig. 1), only a quarter-panel was modeled. Figure 5 shows the quarter-panel finite- 
element model constructed for the sandwich panel. The SPAR constraint commands, SYMMETRY 
PLANE = 1 and SYMMETRY PLANE = 2, were then used to generate the whole panel for thermostruc- 
tural analysis. The panel face sheets were modeled with E43 elements (quadrilateral membrane and bend- 
ing elements), and the sandwich core was modeled with a single layer of S81 elements (hexahedron or 
brick elements) that connect to the upper- and lower-face-sheet E43 elements. 

For the 4S fixed edge (fig. 6(a)) and 4S free edge conditions, the four edges must rotate freely with 
respect to the corresponding edges of the middle plane. To simulate the 4S boundary condition, pin- 
ended rigid rods were attached to the panel edge to connect the two face sheets. The midpoints of these 
rigid rods were pin-jointed to points (fixed or movable in the x- and y-directions) lying in the hypothetical 
middle plane (fig. 6(a)). Each pin-ended rigid rod was modeled with two identical E22 elements (beam 
element for which the intrinsic stiffness matrix is given). To simulate the rigidity of the rods, extensional 
and transverse shear stiffnesses of the E22 elements were made very large. The pin-jointed condition at 
the face sheet edges was simulated by assigning zero values to the rotational spring constants in the stiff- 
ness matrix for the E22 elements. The pin-jointed condition at the hypothetical middle-plane points was 
simulated by eliminating the three rotational constraints. One node of each E22 element was connected to 
the associated node of E43 element, and the other node was connected to the hypothetical middle-plane 
point. The quarter-panel model (fig. 5) for the 4S fixed- and free-edge conditions (to be called 4S model) 
has 1,299 joint locations (JLOCs), 98 E22 elements for the edge rigid rods, 1,152 E43 elements for the 
face sheets, and 576 S81 elements for the sandwich core, as shown in the figure. 

For the 4C fixed-edge condition (fig. 6(b)), the E22 elements at the panel edges may be neglected. 
However, when the E22 elements were attached at the panel edges (i.e., using the 4S model) and enforced 
the zero-edge slopes, the finite-element solutions remained the same as those in which the E22 elements 
were not used. Retaining the E22 elements requires added computational penalty. For the 4C free-edge 
condition, the E22 elements were attached at the panel edges to enforce zero-edge slopes and allow free 
in-plane translations. 


Numerical Input Values 

The dimensions and the material properties used in this study (tables 1 through 3) are identical to the 
titanium honeycomb-core sandwich panel previously used in thermostructural simulation tests at NASA 
Dryden Flight Research Center (ref. 9). 
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Table 1. Geometry of a panel. 


a = b = 24 in 
h = 0.75 in. 
t s = 0.06 in. 


Table 2. Face sheet properties. 



200 °F 

900 °F 

E, lb/in 2 

15.4 x 10 6 

13.1 x 10 6 

V 

0.31 

0.31 

a, in/in-°F 

4.3 x 1CT 6 

5.35 X lcr 6 


Table 3. Honeycomb core (properties at 600 °F). 

E cx = 2.7778 xlO 4 lb/in 2 

E cy = 2.7778 x 10 4 lb/in 2 

E cz = 2.7778 xlO 5 lb/in 2 

G cxy = 0.00613 lb/in 2 

G cy Z = 0.81967 xlO 5 lb/in 2 

G CJ . Z = 1.81 xlO 5 lb/in 2 

v c ^ = 0.658 x Hr 2 

v cyz = 0.643 xlO" 6 

v cxz = 0.643 xlO" 6 

a Q = 5.37 x 10 -6 in/in-°F 

a = 5.37 x 10 s in/in-°F 

a cz = 5.37 xlO -6 in/in-°F 

p hc = 3.674 x 10 -3 lb/in 3 


The temperature loadings on the sandwich panel are T = 900 °F for the upper face sheets (entire 
regions (fig. 2) or central regions (fig. 3)) and T l = 200 °F for the lower face sheet. 

DISPLACEMENT FIELD 

A theoretical equation describing the displacement field of a solid rectangular plate was modified to 
make it applicable to sandwich panels. This equation was required to evaluate the transverse shear effect 
on the sandwich panel deflection. 

The following equation, taken from reference 11, describes the deflection field of a simply supported 
solid isotropic rectangular plate under differential heating. 
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where the origin of the coordinates for this equation is at x = -a/2, y = 0; namely. 


x' = x + - 
2 


( 1 ) 


( 2 ) 


When the transverse shear effect of the sandwich core is neglected (i.e., y = y = 0), the sandwich 
panel behaves like a solid plate, and therefore, the above equation could be usecl to approximate the 
deflection field of the sandwich plate. The validity of equation (1) is addressed later in the “Results” sec- 
tion. To apply .equation (1) for the sandwich plate, the thermal bending term a AT must be replaced with 


a AT = a u T u -a l T l 


( 3 ) 


The displacement field calculated from equation (1) using equation (3), is then compared with that 
calculated from the finite-element method for only the case when the transverse shear effect of the sand- 
wich core is neglected. 


RESULTS 

Displacements 

Figures 7(a) and 7(b) show the half-panel plots of the deformed shapes of the sandwich panel under 
different edge conditions subjected to flat temperature profile heating. These half-panel plots were gener- 
ated from the quarter-panel plots by using the SYMMETRY command. The panel’s deformed shapes 
under fixed and free edges are the same; however, as will be seen later, the induced thermal stress fields 
are quite different. Notice that under the 4C edge condition (fixed or free, fig. 7(b)), the panel deflection 
is zero (i.e., the deformed panel remains flat). 

Figures 8(a) and 8(b) show the similar half-panel plots for the case when the heating is of the dome- 
shaped temperature profile. Again, the deformed shapes for the fixed- and free-edge cases are identical. 
Unlike the previous case, the panel deflection away from the boundaries under the 4C edge condition, 
fixed or free, (fig. 8(b)) is nonzero. 

Figure 9 shows the deflection curves of the sandwich panel’s middle plane center line, along the 
x-axis, for different heating cases. The figure also shows the deflection curve calculated from 
equation (1) up to 10 terms summation for flat temperature profile heating. The deflection curve calculat- 
ed from equation (1) falls pictorially on that calculated from the finite-element method (for the flat tem- 
perature case) neglecting the transverse shear effect (i.e., by setting y xz = y yz = 0, namely by making 
G* and G„ v _ very large). The maximum deflection at the center of sandwich panel middle plane w 

CA,Z L-y i r max 

calculated from equation (1) using 10-terms series summation is w = 0.2931 in., and that calculated 
from the finite-element method (for y xz = y yz = 0) is w max = 0.2920 in. The close correlation between 
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these two values gives confidence in the applicability of equation (1) for the sandwich panel, and also in 
the adequacy of the finite-element model used. 

The deflection equation (1) and the finite-element models were also used to study the effect of trans- 
verse shear on panel deflections. Table 4 lists the panel maximum deflections w max for different heating 
and edge conditions. The w max values in parentheses are applicable when the transverse shear effect 
is neglected (i.e., y xz = y yz - 0). The w max value calculated from equation (1) (10 terms summation) is 
also shown in brackets for comparison. Notice that by neglecting the transverse shear effect, the panel 
deflection could be underpredicted by 3 to about 11 percent depending on the edge and heating condi- 
tions. In practical application, the panel is under dome temperature profile heating and 4C edge condition 
(closer to fixed edges rather than free edges). For this case, the maximum deflection is only w max = 
0.0576 in. Such a small panel deflection greatly minimizes any concern that the deformed panel could 
severely disturb the airflow field and, therefore, alter the surface heating rate. 


Table 4. Maximum deflections at center of sandwich panel’s middle plane; 
T u = 900 °F, T t = 200 °F. 




w 

max 7 

in. 


Temperature profile 

4S fixed 

4S free 

4C fixed 

4C free 


0.3305 

0.3305 

0 

0 

Flat 

(0.2920) 

(0.2920) 

— 

— 


[0.2931] 

— 

— 

— 

Dome 

0.3050 

0.3050 

0.0576 

0.0576 

(0.2745) 

(0.2745) 

(0.0557) 

(0.0557) 


Thermal Stresses 

This section presents thermal stress results for flat and dome-shaped temperature heating applied to 
the sandwich panel under the various edge conditions. 

Flat Temperature Profile 

Upper Face Sheet 

Figures 10 through 21 show various distributions of normal stresses {a x , o y }, which are negative 
(i.e., compression), and the shear stress (x ), in the upper face sheet of the sandwich panel induced by 
flat-temperature-profile heating. Again, these half-panel plots were generated from quarter-panel plots 
using the SYMMETRY command. The figures show the peak stress points and values of peak stress. For 
the 4S fixed- and free-edge cases, the distributions of {o x , o^} (figs. 10, 11, 13, and 14) are slightly sad- 
dle shaped, and the peak compression points of { o x , o y } are at the midpoints of the panel edges, y = b/2 
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and x = a/2, respectively. The distributions of the shear stress x in the upper face sheet for the 4S fixed- 
and free-edge cases (figs. 12 and 15) are distorted bell-shaped within each quarter-panel region. The mag- 
nitude of x reaches its peak value near the panel corners and decreases steeply to zero at the panel edges 
and rapidly to zero at the two axes of symmetry. By freeing the panel edges from the fixed constraint, the 
magnitudes of compressive stresses {o x , o y } could be reduced considerably (figs. 10, 11, 13, and 14); 
however, the distributions of shear stress x xy remained almost the same (figs. 12 and 15). For the 4C 
fixed-edge case, the compressive stresses {c ^.,'0 } in the upper face sheet (figs. 16 and 17) are constant 
everywhere, and the shear stress x is zero everywhere (fig. 18). For the 4C free-edge case, the compres- 
sive stresses {o x , o y } are almost constant over the upper face sheet (figs. 19 and 20), and the shear stress 
x induced in the upper face sheet (fig. 21) is nearly zero. 

Lower Face Sheet 

Figures 22 through 33 show various distributions of {o x , a , x } induced in the lower face sheet of 
the sandwich panel under different edge conditions for flat- temperature-profile heating. For the fixed- 
edge cases (4S and 4C), {o x , o y } are negative (i.e., compression); however, for the free-edge cases (4S 
and 4C), {o^, o y } are positive (i.e., tension). For the 4S fixed-edge case, the peak compression points of 
{o^, o y } (figs. 22 and 23) are no longer located at the edge midpoints like the upper-face- sheet case 
(figs. 10 and 11). For the 4S free-edge case, similar to the upper-face-sheet case for which {a , a } are 
negative (figs. 13 and 14), the peak tensile stress points of {o^, a } (figs. 25 and 26), are at the midpoints 
of the panel edges. The distributions of shear stress x yv in the lower face sheet for the 4S fixed- and free- 
edge cases (figs. 24 and 27) are very similar to those for the upper-face-sheet case (figs. 12 and 15); 
however, the sign of x is reversed. For the 4C fixed-edge case, the compressive stresses {a^, a^} are 
identical and constant everywhere in the lower face sheet (figs. 28 and 29), but the magnitude is much 
lower than that in the upper face sheet (figs. 16 and 17), and the shear stress x in the lower face sheet 
(fig. 30), like the upper face sheet case (fig. 18), is zero everywhere. For the 4C free-edge case, {a x , o y } 
in the lower face sheet are both positive and nearly constant (figs. 31 and 32), and the shear stress x 
induced in the lower face sheet (fig. 33) is at an insignificant level similar to the case for the upper face 
sheet (fig. 21). 

Sandwich Core 

Figures 34 and 35 show the distributions of transverse shear stresses { x xz , x } induced in the sand- 
wich core under 4S fixed- or 4S free-edge condition subjected to flat-temperature-profile heating. Each 
transverse shear stress value used in the plots is the average of the eight stress values at the eight nodes of 
the S81 element. Both the 4S fixed- and 4S free-edge cases induced identical transverse shear stresses. 
The shape of x xz plot (fig. 34) is like an airplane wing with winglets; for x yz (fig. 35) the shape is like 
fox ears near the panel’s edge. The values of {x xz , x yz } are quite low in the core central region and rises 
steeply to their respective peak values at the comers of the core (i.e., the transverse shear stress concen- 
trations occur at the panel’s comers). The levels of {x xz , x yz } stress concentrations are relatively low; 
however, they might cause the thin honeycomb walls to buckle in shear. 

Under flat-temperature-profile heating, both the 4C fixed- and 4C free-edge cases induced no 
transverse shear stresses {x xz , x yz } in the sandwich core, as shown in figures 36 and 37. For actual 
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application, the panel’s edge condition is closer to the 4C fixed-edge condition, which induces no trans- 
verse shears. Therefore, under the present heating level, no concern about stress concentration caused by 
transverse shear is warranted. 

Dome Temperature Profile 

Upper Face Sheet 

Figures 38 through 49 show the distributions of {a^, cy x } in the upper face sheet of the sandwich 
panel induced by the dome-temperature-profile heating. The distributions of the compressive stresses 
{o x , o y } for the 4S fixed- and free-edge cases (figs. 38, 39, 41, and 42) and 4C fixed- and free-edge cases 
(figs. 44, 45, 47, and 48), are hat shaped, reflecting partly the shape of the dome temperature profile. 
Unlike the flat-temperature-profile case, the peak compression points of { o x , o y } are now at the bound- 
ary of the central plateau zone and not at the panel’s edges. The shear stress distributions for the 4S fixed- 
and free-edge cases (figs. 40 and 43) are similar to those of the flat-temperature-profile case (figs. 12 and 
15) but with slightly higher peak magnitudes. For the 4C fixed- and free-edge cases, the shear stress x xy 
(figs. 46 and 49), is zero only at the two axes of symmetry and the panel’s comers (the plotted comer 
points are not exactly at the panel’s comers, and therefore, show finite values) but nonzero at the panel 
edges because of the nonuniform thermal expansions. The peak magnitudes of x are lower than those 
for the 4S fixed- and free-edge cases (figs. 40 and 43) and are in the vicinity of the panel’s corners. 

Lower Face Sheet 

Figures 50 through 61 show various distributions of { o x , a , x } in the lower face sheet induced by 
the dome-temperature-profile heating. The stress distributions are very similar to those for the flat tem- 
perature profile and, unlike the upper face sheet, they do not reflect the dome temperature profile. Those 
figures also indicate the peak stress points. For the 4C fixed and 4C free cases (figs. 56 through 61), the 
distributions of {a , o , x } stresses in the lower face sheets are slightly wavy in shape with peak stress 
points at the edges. 

Sandwich Core 

Figures 62 and 63, respectively, show the distributions of transverse shear stresses {x xz , x yz } induced 
in the sandwich core for the 4S fixed- and 4S free-edge cases under dome-temperature-profile heating. 
Both the 4S fixed- and 4S free-edge case give identical { x xz , x yz } distributions. The x xz plot (fig. 62) 
also looks like an airplane wing with winglets, and the x plot (fig. 63) looks like cat ears near the 
panel’s edge. The peak magnitudes of { x xz , x yz } are now at the core edges and near the core comers — not 
exactly at the core corners like the previous case. The {x xz , x yz } stress concentrations are less severe for 
the dome-temperature-profile case (figs. 62 and 63), as compared with the flat- temperature-profile case 
(figs. 34 and 35). For dome-temperature-profile heating, the magnitudes of transverse stresses {x xz , x yz } 
induced in the sandwich core under both 4C fixed- and 4C free-edge conditions (figs. 64 through 67) are 
quite close and very small — not exactly zero like for flat-temperature-profile heating (figs. 36 and 37). 
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Peak Stress Summary 

Table 5 lists the peak values of the thermal stresses { o x , o , x } (positive or negative) induced in the 
upper and lower face sheets of the sandwich panel. Table 6 lists the peak values of the transverse shear 
stresses {x X2 , x yz j induced in the sandwich core. 

Table 5. Peak thermal stresses in face sheets of sandwich panel; T = 900 °F, T, = 200 °F (yield stress = 
126,000 lb/in 2 ). “ 


Upper face sheet 


Edge 

condition 

Flat temperature profile 

Dome temperature profile 

o x , lb/in 2 

o y , lb/in 2 

x xy , lb/in 2 

o x , lb/in 2 

o y , lb/in 2 

"V lb/in2 

4S fixed 

-89,245 

-87,430 

23,827 

-75,961 

-75,796 

28,579 

4S free 

-26,161 

-25,783 

23,904 

-26,140 

-25,576 

30,917 

4C fixed 

-107,470 

-107,470 

0 

-96,192 

-96,087 

12,984 

4C free 

-46,177 

-45,801 

219 

-45,488 

-44,996 

14,062 

Lower face sheet 

4S fixed 

-62,948 

-64,010 

23,827 

-53,433 

-54,553 

19,242 

4S free 

25,671 

22,931 

23,755 

21,084 

18,251 

17,086 

4C fixed 

-19,190 

-19,190 

0 

-29,399 

-28,160 

2,115 

4C free 

44,339 

44,393 

219 

36,994 

36,743 

2,757 


Table 6. Peak transverse shear stresses in sandwich core; T u = 900 °F, T t = 200 °F. 


Edge 

condition 

Flat temperature profile 

Dome temperature profile 

x xz , lb/in 2 

x yz , lb/in 2 

Ib /> n2 

x yz , lb/in 2 

4S fixed 

2,216 

1,798 



4S free 

2,216 

1,798 



4C fixed 

0 

0 

134 

86 

4C free 

0 

0 

142 

88 


CONCLUSIONS 


Finite-element thermal stress analyses were performed on a titanium honeycomb-core sandwich 
panel supported at its edges under four different edge conditions, and heated on one side under both flat 
and dome-shaped temperature profiles. Detailed deformation and thermal stress fields induced in the 
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sandwich panel were presented graphically for easy visualization. The key results of the thermostructural 
analyses are as follows: 

1. If the transverse shear effect of a sandwich core is neglected, the maximum deflection of the 
sandwich panel was underpredicted by 3 to about 11 percent, depending on edge conditions and 
heating temperature profiles. 

2. Under the flat-temperature-profile heating: 

(a) The classical deflection equation for a simply supported rectangular flat plate adequately 
describes the deformation field of a simply supported sandwich panel for which the transverse 
shear effect of the sandwich core is neglected, as validated by the finite-element solutions. 

(b) For the 4S fixed- and 4S free-edge conditions, the peak stress points of the normal stresses 

{o x , o } are at the edges of the face sheets; the peak stress points of the shear stress x are at 
y xy 

the diagonal lines and near the comers of the face sheets; and the peak stress points of the 
transverse shear stresses {x xz , x yz } are right at the comers of the sandwich core. 

(c) For the 4C fixed-edge condition, the normal stresses (a , a } in both of the face sheets are 

x y 1 

constant everywhere, and the shear stress x xy in both of the face sheets is zero everywhere. 

For the 4C free-edge condition, the normal stresses {a , or } in both of the face sheets are 

x y 

almost constant, and the shear stress x xy in both of the face sheets is negligibly low. For both 
4C fixed- and 4C free-edge conditions, the transverse shear stresses {x xz , x z } are zero every- 
where in the sandwich core. 

3. Under the dome-temperature-profile heating: 

(a) For all four edge conditions, the peak stress points of normal stresses {<j x , o y } in the upper 
face sheet are at the boundary of the temperature plateau zone — not at the face sheet edges. 
The peak stress points of the shear stress x in the upper face sheet (for all the edge condi- 
tions) and in the lower face sheet (for 4S fixed- and 4S free-cases only) are on the diagonal 
lines and near the comers of the face sheets. 

(b) The distributions of {a x , o y } in the lower face sheet for the 4S fixed and 4S free cases are 
very similar to those for the flat-temperature-profile case and do not reflect the temperature 
profile applied to the upper face sheet. For the 4C fixed and 4C free cases, the distributions of 
{o x , o y , x } stresses in the lower face sheets are slightly wavy in shape with peak stress 
points at the edges of the lower face sheet. 

(c) For 4S fixed and 4S free cases, the peak stress points of the transverse shear stresses {x xz , 
v yz } in the sandwich core are at the core edges and very near the core comers. For the 4C 
fixed and 4C free cases, the transverse shear stresses {* xz , * yz } induced in the sandwich core 
have extremely low magnitudes. 


Dry den Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, May 3, 1996 
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Figure 4. Measured temperature distribution in upper surface of titanium honeycomb-core sandwich panel, 
heated on upper side at 10 °F/sec heating rate, with four edges supported by test fixtures (heat sink). 
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(a) 4S edge condition (fixed and free). 
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(b) 4C edge condition (fixed and free). 

Figure 8. Deformed shapes of sandwich panel under heating on one side; dome temperature profile; T 
900 °F, T t = 200 °F; half-panel plots. 
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Figure 9. Deflections of sandwich panel’s middle plane along x-axis; T u = 900 °F, T ^ = 200 °F. 
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Normal and Shear Stress Distributions in Upper Face Sheet — Flat Temperature Profile 
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Figure 10. Distribution of o x in the upper face sheet; 4S fixed-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200. °F. 
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Figure 11. Distribution of a in the upper face sheet; 4S fixed-edge condition; flat temperature profile; 
T = 900 °F, T, = 200 °F. V 
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Figure 12. Distribution of % x in the upper face sheet; 4S fixed-edge condition; flat temperature profile; 
T = 900 °F, T, = 200 °F. 

U * 


22 






2 


in 


x = 23,904 lb/in 2 


Panel 

corner 











T 


T 


[illil 


Panel 


corner 



Figure 16. Distribution of a in the upper face sheet; 4C fixed-edge condition; flat temperature profile; 


T u = 900 °F; T t = 200 °F. 
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Figure 17. Distribution of a in the upper face sheet; 4C fixed-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 18. Distribution of % x in the upper face sheet; 4C fixed-edge condition; flat temperature profile; 
T.. = 900 °F. T, = 200 °F. ^ 
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Figure 19. Distribution of o x in the upper face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 20. Distribution of a in the upper face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T { = 200 °F. 
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Figure 21. Distribution of x in the upper face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T { = 200 °F. 
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Normal and Shear Stress Distributions in Lower Face Sheet — Flat Temperature Profile 
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Figure 28. Distribution of o x in the lower face sheet; 4C fixed-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 29. Distribution of o in the lower face sheet; 4C fixed-edge condition; flat temperature profile; 
T u = 900 °F, T x = 200 °F. 
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Figure 30. Distribution of x in the lower face sheet; 4C fixed-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 31. Distribution of o x in the lower face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 



Figure 32. Distribution of a in the lower face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 33. Distribution of % in the lower face sheet; 4C free-edge condition; flat temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Transverse Sheer Stress Distributions in Sandwich Core — Flat Temperature Profile 
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Figure 34. Distribution of transverse shear stress x xz in sandwich core; 4S fixed- or 4S free-edge condi- 
tion; fiat temperature profile; T u - 900 °F, T l - 200 °F. 
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Figure 35. Distribution of transverse shear stress x in sandwich core; 4S fixed- or 4S free-edge condi- 
tion; flat temperature profile; T u = 900 °F, T l = 200 °F. 
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Figure 36. Distribution of transverse shear stress x in sandwich core; 4C fixed- or 4C free-edge condi- 
tion; flat temperature profile; T u = 900 °F, T l = 200 °F. 



Figure 37. Distribution of transverse shear stress x yz in sandwich core; 4C fixed- or 4C free-edge condi- 
tion; flat temperature profile; T u = 900 °F, T [ = 200 °F. 
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Normal and Shear Stress Distributions in Upper Face Sheet — Dome Temperature Profile 
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Figure 38. Distribution of o x in the upper face sheet; 4S fixed-edge condition; dome temperature profile; 
T u = 900 °F, T ( = 200 °F. 
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Figure 39. Distribution of a in the upper face sheet; 4S fixed-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. y " 
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Figure 40. Distribution of x in the upper face sheet; 4S fixed-edge condition; dome temperature 
profile; T u = 900 °F, T { = 200 °F. 
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Figure 44. Distribution of a in the upper face sheet; 4C fixed-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 45. Distribution of o in the upper face sheet; 4C fixed-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 46. Distribution of % in the upper face sheet; 4C fixed-edge condition; dome temperature 
profile; T u = 900 °F, = 200 °F. 
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T u = 900 °°F, T l = 200 °F. 
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Normal and Shear Stress Distributions in Lower Face Sheet — Dome Temperature Profile 
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Figure 50. Distribution of cs x in the lower face sheet; 4S fixed-edge condition; dome temperature profile; 
T u = 900 °F, T { = 200 °F. 
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Figure 5 1 . Distribution of o in the lower face sheet; 4S fixed-edge condition; dome temperature profile; 
T u = 900 °F, T { - 200 °F. V 
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Figure 52. Distribution of x in the lower face sheet; 4S fixed-edge condition; dome temperature 
profile; T u = 900 °F, T { = 
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Figure 53. Distribution of o x in the lower face sheet; 4S free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 54. Distribution of a in the lower face sheet; 4S free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 55. Distribution of x in the lower face sheet; 4S free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 56. Distribution of o x in the lower face sheet; 4C fixed-edge condition; dome temperature profile; 
T u = 900 °F; T t = 200 °F. 
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Figure 57. Distribution of a in the lower face sheet; 4C fixed-edge condition; dome temperature profile; 
T u = 900 °F, T { = 200 °F. y 
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Figure 58. Distribution of x in the lower face sheet; 4C fixed-edge condition; dome temperature 
profile; T u = 900 °F, T t = 200 °F. 
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Figure 59. Distribution of o x in the lower face sheet; 4C free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 60. Distribution of o in the lower face sheet; 4C free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Figure 61. Distribution of x in the lower face sheet; 4C free-edge condition; dome temperature profile; 
T u = 900 °F, T t = 200 °F. 
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Transverse Shear Stress Distributions in Sandwich Core — Dome Temperature Profile 
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Figure 62. Distribution of transverse shear stress x xz in sandwich core; 4S fixed- or 4S free-edge condi- 
tion; dome temperature profile; T u = 900 °F, T l = 200 °F. 




Panel 

corner 


960503 

Figure 63. Distribution of transverse shear stress x yz in sandwich core; 4S fixed- or 4S free-edge condi- 
tion; dome temperature profile; T u = 900 °F, T [ - 200 °F. 
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Figure 64. Distribution of transverse shear stress t xz in sandwich core; 4C fixed-edge condition; dome 
temperature profile; T u = 900 °F, T l = 200 °F. 



960505 

Figure 65. Distribution of transverse shear stress x in sandwich core; 4C fixed-edge condition; dome 
temperature profile; T u = 900 °F, T l = 200 °F. 
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Figure 66. Distribution of transverse shear stress t xz in sandwich core; 4C free-edge condition; dome 
temperature profile; T u - 900 °F, T t = 200 °F. 
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Figure 67. Distribution of transverse shear stress % in sandwich core; 4C free-edge condition; dome 
temperature profile; T u - 900 °F, T t = 200 °F. 
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